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Abstract. The peculiar dip in the outer rotation curve at a distance of 9 kpc, which 
was recently confirmed by precise measurements with the VERA VLBI array in Japan, 
suggests donut-like substructures in the dark matter (DM) halo, since spherical or elliptical 
distributions will not cause a dip. Additionally, such a donut-like DM structure seems to be 
required by the dip in the gas flaring of the disk. In this paper we consider the impact of such 
DM substructure in the disk on the rotation curve, the gas flaring, the local DM density and 
the local surface density. A global fit shows that the rotation curve is best described by an 
NFW DM profile complemented by two donut-like DM substructures at radii of 4.2 and 12.4 
kpc, which coincide with the local dust ring and the Monocerus ring of stars, respectively. 
Both regions have been suggested as regions with tidal streams from "shredded" satellites, 
thus enhancing the plausibility for additional DM. If real, the radial extensions of these 
nearby ringlike structures enhance the local dark matter density by a factor of four to about 
1.3±0.3 GeV/cm 3 . We find that i) this higher DM density is perfectly consistent with the 
local gravitational potential determining the surface density and ii) the s-shaped gas flaring 
is explained. Such a possible enhancement of the local DM density is of great interest for 
direct DM searches and the ringlike structure would change the directional dependence of 
gamma rays for indirect DM searches. 
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1 Introduction 

The flat rotation curves (RC) of galaxies show that dark matter (DM) is present in the form 
of large scale halos with a density slowly falling off from the center (like 1/r 2 for a flat RC), 
see e.g. [1, 2]. In our Galaxy the rotation curve is not flat but has two peculiar dips, one at 3 
kpc and one at 9 kpc [3]. These dips can not be explained by smoothly varying components, 
but need ringlike perturbations in the disk, as studied in detail in Ref. [3]. These authors 
combined all available rotation curve data and assumed substructure in the visible sector to 
explain the dips. However, especially for the "outer" ring located outside the solar radius 
of 8.3 kpc, the amount of visible matter rapidly decreases and it seems hard to explain the 
increase in RC with a sudden increase in visible matter. 

To determine the local DM density one has usually ignored this increase in the RC, 
assuming that the distance measurements are too uncertain, see e.g. [4]. The results for the 
local DM density are then around 0.3-0.4 GeV/cm 3 [5-7] with the precise value depending 
on the assumptions and fitting procedures. E.g. in Ref. [6] it is shown that the various disk 
parameters are almost 100% correlated with the DM profile parameters, which might be at the 
origin of the rather large differences in error estimates. A local DM density of 0.3 GeV/cm 3 
is only 8% of the local matter density, as determined from the precise measurements of the 
gravitational potential using the star counts and velocity measurements from the Hipparcos 
satellite [8-10] in agreement with earlier measurements of the surface density [11-13]. Most 
of these authors suggested from the decrease in the inner RC combined with the low surface 
density that there is no need for a significant contribution of DM in our solar neighborhood. 
However, there are four additional new pieces of evidence, that bear on this question: 

• Measurements with the VERA VLBI array in Japan using the maser light from molec- 
ular clouds confirm the peculiar change of slope around 9 kpc [14]. With the precise 
distance measurements one cannot ignore the high rotation speed in the outer RC 
anymore. 

• The thickness of the gas layer increases with increasing distance from the Galactic 
center, the so-called gas flaring. The thickness is given by the gravitational potential 
and shows a clear "dip" in the outer Galaxy, which seems only explainable by a donut- 
like DM substructure in the Galactic plane [15]. 
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• The Monocerus ring of stars in the outer Galaxy was originally discovered in the SDSS 
data [16, 17]. Subsequently it was found to span an almost circular structure of low- 
metallicity stars of at least 170 degrees around the Galaxy in the plane of the disk, see 
e.g. [18-22]. Such a large overdensity over a large range is suggestive of the infall of a 
dwarf galaxy and N-body simulations using the Monocerus ring as input predict DM 
ringlike structures in the Galactic plane with the closest and most dense ring located 
at 13 kpc from the Galactic center, see Refs. [23-25] and references therein. This leads 
to DM densities in the plane of the disk in broad agreement with the gas flaring, which 
forms a strong argument against interpretations of the Monoceros ring as extensions of 
the disk, which led to lively discussions, see e.g. [26-28]. 

• additional dark matter in the inner region has been suggested from the early infall of 
a dwarf galaxy, presumably some 13 billion years ago [29-34] , which created the thick 
disk. The thick disk has a significantly larger scale height and its mainly metal-poor 
stars imply a different population of old stars. A possible hint of a ringlike gravitational 
potential well in the inner Galaxy is provided by the ringlike structure of the dust ring 
with the highest density of the dust ring located at a radius of about 4 kpc from the 
Galactic center and a radial extension of about 1 kpc [35] . 

It is the purpose of this paper to reinvestigate the local DM density by assuming two 
donut-like overdensities of matter in the Galactic plane, one inner ring and one outer ring 
in order to explain the dips in the rotation curve at 3 and 9 kpc as well as the dip in the 
gas flaring. The data on the RC from Ref. [3] was supplemented with the new precise 
measurements with the VERA VLBI array in Japan using the maser light from molecular 
clouds [14]. The fitting procedure is described in our previous paper, where the new VERA 
data were not available yet and the dip in the rotation curve at 9 kpc was ignored [6]; this 
paper will be referred to as Paper I. 



2 Analysis 

2.1 Rotation curve 

The DM profile including the rings can be parameterized as follows: 



7-/3 




(2.1) 



where the first term in Eq. 2.1 represents the smooth Galactic halo and the second term 
describes the DM donut-like rings. For the smooth DM halo an NFW profile with a = 1, (5 = 
3, 7 = 1 in Eq. 2.1 and a scale radius of 12 kpc was taken, which yields an excellent fit to the 
inner rotation curve [6] , as displayed by the short-dashed line in Fig. 1 . In the rings the DM 
distribution in z-direction is taken to decrease exponentially with the scale height a z n and 
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Figure 1. The Galactic rotation curve (RC). The dashed line corresponds to a fit of the inner RC 
only with a smooth NFW profile. The solid line includes DM ring-like substructure at 4.2 and 12.4 
kpc. For clarity the weighted average of many data points has been plotted, where the data was taken 
from Ref. [6] with in addition the precisely measured ones with VBLI technology from Ref. [14], 
shown as upside-down (red) triangles. 



distributed in r-direction around the ring radius R n with a Gaussian width a^ n . The rings 
have been assumed to be circular, i.e. the ellipticities e xy and e z have been set to 1, although 
in reality infall of dwarf galaxies will never lead to perfectly circular rings. However, the 
used data is most sensitive to the nearest segment of the stream, which seems locally well 
described by a circular structure and changing the ellipticity did not improve the fit. The 
Gaussian radial width of the outer ring is taken to be slightly different for the inner and the 
outer side 1 , which is required for the steep increase of the slope in the rotation curve after 
9 kpc combined with the slow decrease at larger radii. Such a difference is expected from 
the conservation of angular momentum for tidal streams, since they cannot reach arbitrary 
small Galactocentric distances. The outer RC is much better described, when the rings 

1 The radial density was taken to decrease at the inner side as an s-shape built from two parabolic functions, 
which reach zero at 5 kpc from the center of the ring, i.e. at R=7.4 kpc. 
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p n [GeV cm~ 3 ] 
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o> [kpc] 


2.5 


3.2 


a z [pc] 
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Table 1. Ring parameters of the best \ 2 fit- The errors are typically 30%, but are correlated with 
the disk parameters, which have been fixed to the values given in Ref. [6]. 
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Figure 2. The radial dependence of the DM density including the ringlike structures at 4.2 and 12.4 
kpc. 



are included, as shown by the solid line in Fig. 1. The parameters of the outer ring are 
dominated by the precise VLBI data points from VERA. The main reason for the decreasing 
inner RC in this scenario is not the lack of DM, but the large negative contribution from the 
outer ring, which creates an outward gravitational force, thus reducing the centrifugal force 
and hence the rotation speed inside the ring. The ring parameters for the best fit to the data 
are given in Table 1. 
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These parameters provide DM rings with masses of 1.4 • 10 solar masses for the inner 
ring and 3.9 • 10 10 solar masses for the outer ring. Although these masses are only a minor 
contribution to the total mass of about 10 12 solar masses, one may worry about the disk 
stability and probability of accreting such heavy masses. A recent paper on the assembly 
history of DM halos [36] finds that in three of their six simulated galaxies the central mass 
was brought in through mergers involving host halos with mass greater than 10 10 /i _1 solar 
masses in agreement with previous studies [25]. Furthermore they state that "late-accreted 
mass typically settles in the outskirts of a halo, thus allowing the hierarchical growth of halos 
to be reconciled with the ubiquitous presence of thin stellar disks" , again in agreement with 
Ref. [25]. Also the thin disk may be regrown from the cooling of the heated gas after a 
merger, see e.g. [37] and if the disk survival depends sensitively on the way how the accreted 
mass approaches the Galaxy and how it is introduced in the simulations [33]. It should be 
remembered that our neighboring galaxy, Andromeda (M31), shows also a ringlike structure 
in the outer galaxy, which is attributed to the infall of a dwarf galaxy as well [38] . 

The local density of the smooth NFW halo P0,Haio is found to be 0.28 GeV/cm 3 , while 
the inner and outer rings contribute 0.87 and 0.14 GeV cm~ 3 , respectively, which leads to: 

p DM = 1.3 ± 0.3 GeV/cm 3 = 0.030 ± 0.007 M /pc 3 , (2.2) 

so the rings enhance the local DM density by about a factor four. The error is systematic 
and depends mainly on the assumed parameterizations of the rings. The radial dependence 
of the DM density within the Galactic disk is shown in Fig. 2. The influence of the DM rings 
is clearly visible. 

2.2 Surface density 

The enhanced DM density from the rings leads to a different gravitational potential in the 
solar neighborhood ("local" potential), which in turn is expected to change the spatial dis- 
tribution and velocity distribution of the local star populations in much the same way as the 
gas distribution in the earth's atmosphere is determined by its gravitational potential. So the 
distribution of star counts as function of galactic height above the disk yields information 
on the total local matter density, which can be compared with the visible matter density, 
thus yielding information on the local DM density. The principle was pioneered by Jan Oort 
in the 1960's, but a seminal study using much better data was performed by Kuijken &; 
Gilmore (KG) [11-13]. They concluded that locally little dark matter was needed, a con- 
clusion confirmed later by independent analyses using the precise data from the Hipparcos 
satellite [9, 10]. These data provided also a precise measurement of the local mass density 
(Oort limit), which includes visible and dark matter [8, 9]: 

p Qjtot (z = 0) = 0.102 ± 0.005 M Q pcT 3 . (2.3) 

From a comparison with Eq. 2.2 one observes that the local DM density is about 30% of the 
total density, which is clearly a significant contribution. How can it be that previous authors 
concluded from the same star counts that little DM is needed? 

To understand this surprising result, let us go back and see what is really measured. 
The number density and velocities of a tracer population is measured as function of the 
distance above the disk. We assume that the stars move on nearly circular orbits in the 
plane, so the vertical motions are independent of the circular motions. Furthermore, we 
assume that the tracer population is well-mixed, i.e. the phase-space distribution f(z,v z ) 
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Figure 3. The vertical gravitational potential at the position of the Sun $(z) for a fit with and 
without DM ringlike substructures. The data points (derived from the data in Ref. [10]) and its 
errors are shown only for comparison, since these data points have not been used in the fit to the 
constraints mentioned in Paper I. The errors are taken to be 30%, as estimated from the uncertainties 
in the phase-space distribution of the Hipparcos data and the parametrization of the potential. 



is in a steady-state in a static potential. Then the energy at the local solar Galactocentric 
radius Rq can be written as E z = &(Ro, z ) + \ v 1 an d is an integral of motion with E z < for 
bound stars. The distribution function is given by the Boltzmann factor: f(z, v z ) = f(E z ) oc 
ex.p(—E z /a(z) 2 ), where the velocity dispersion is defined by <r(z) 2 =< v 2 > — < v z > 2 . 
Integrating over the velocity yields for the number density of stars as function of z: n(z) = 
no exp(— &(Ro, z)/a(z) 2 ). So measuring the number density n(z) and the velocity dispersion 
<r(z) of stars as function of z determines the local gravitational potential <&(Ro, z) as function 
of z. 

To obtain the potential from the data it is convenient to have an analytical parametriza- 
tion. KG modeled the potential as the sum of two components [11]: 

$ z /2ttG = S (\/z 2 + D 2 - £)) + Pc{i z 2 . (2.4) 

The first term is the potential of an exponential disk with a scale height D and density So, 
while the latter represents the potential of a locally constant halo DM density p c Q. If the 
radial and z movements are independent, p e g equals the local DM density Po,tot(z = 0), which 
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is furthermore assumed to be independent of z for the small volume considered. For precise 
enough data one could obtain the three parameters So, p e ff and D from the distribution 
function. However, these parameters are strongly correlated and the data are not precise 
enough. Therefore additional measurements are welcome. The halo density p e g can be 
obtained from the rotation curve, as discussed above. Furthermore, the Oort limit (Eq. 2.3) 
is given by the sum of the DM density and the disk density: 

ftotal(^ = 0)=p e ff + S /(2 J D). (2.5) 

This yields a relation between So and D, so if p e g is obtained from the RC, one basically has 
only one free parameter left, namely D, which can be precisely determined from a fit to the 
distribution function. Unfortunately, a different choice of p e Q will yield a different value of 
the scale parameter D via the constraint in Eq. 2.5. What basically happens: p{z) is fixed 
at z = by the Oort limit and for z » D by the p e g from the RC and the interpolation at 
intermediate z is given by D. If we go from a RC without rings to a RC curve with rings, 
the values of p e g and D will change. However, the potential in Eq. 2.4, as probed by the 
distribution of stars, may change less and be consistent with the data. To see this one can 
relate an axisymmetric mass distribution to the potential via the Poisson equation: 



<9 2 <3?fi? z) Id 



(2.6) 



The surface density is obtained by integrating the density along z, which is then proportional 
to the first derivative of the potential or the vertical force F z . This force can be written in 
terms of the parametrization of the potential as: 

F z /2irG=-^= + 2p eB . (2.7) 
\J z z + D z 

From the Poisson equation it is clear that the surface density is obtained by integrating the 
potential, so the surface density is proportional to the vertical force plus a small correction 
term, if the force varies with radius: 

2 

tW«W)-/^^-e(".-is)' < 2 - 8 > 

— 2 

Here v r = i?<9<3?(i?o, z)/dR is the local rotation speed. If we choose p e g from the rotation 
curve without rings, we obtain results similar to previous measurements for \z\ = 1.1 kpc: 
Siot (< | z |) = 71 ± 6 Mq pc -2 . However, if we choose the four times higher value of p e g from 
the RC including the rings, we find a considerably higher value: 

S to t (< |*|) = 92 ±7 M &P c~ 2 . (2.9) 

Thus an increase in p e Q leads to a larger surface density, but this changes only the slope of 
the potential for z >> D, not so much its value, as demonstrated in Fig. 3. Here we simply 
plotted the potential as calculated for the mass distribution obtained from a fit to the RC 
with and without rings. For the fit the surface density was not set to the value used in Paper 
I, but became a function of p c s~. S tot = a + b • /O 0) dm with a = 57.5 M Q pc~ 2 and b = 1.1 • 10 3 
pc. This relation was obtained from a fit to the values in Table 3 of Ref. [10]. 
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Figure 4. Gas flaring of a DM halo profile combined with a substructure of two DM rings at 4.5 and 
12.5 kpc. 



To check that we have a self-consistent solution the surface density was calculated first 
by direct integration of the density using the mass model in Paper I for the visible matter 
and Eq. 2.1 for the DM density, i.e. we use the first equality in Eq. 2.8. Secondly, the local 
gravitational potential was calculated from the mass distributions for visble and DM and 
the vertical force was calculated from the derivative of the potential leading to the surface 
density using the last equality in Eq. 2.8. Both methods agree within errors. 

In summary, the distribution function of local stars probe the potential, not the surface 
density. The latter is determined by the slope of the potential and the DM dominates only 
in the region 0.5 < \z\ < 1.1 kpc. A fourfold increase in the DM density changes the slope of 
the potential in this region by about 20%, which is inside the experimental errors. 

2.3 Gas flaring 

From the gravitational potential one can calculate the corresponding gas flaring, which is the 
effect of the increase in the half-width-at-half-maximum (HWHM) of the gas layer as function 
of radial distance from the Galactic center due to the decrease in gravitational potential. 
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Of course, the gas flaring will be affected by the donut-like DM rings. In analogy to the 
barometric equation the vertical decrease of the interstellar hydrogen can be parameterized 
as: 



The experimental determination of the velocity dispersion u yielded values of 7 ± 1 km s 
[39] and 8 ± 1 km s~ x [40]. Fig. 4 shows the gas flaring for a potential with and without DM 
rings, as obtained from the measurement of the spatial distribution of the 21 cm emission 
line of atomic hydrogen. Note that the gas flaring is independent of the warping of the disk, 
since the center-of- gravity of the gas layer was allowed to differ from the galactic plane. The 
gas flaring is different for the northern and southern hemisphere [15], which indicates that 
either the gravitational potential is not axisymmetric in the disk or local phenomena, like 
recent star bursts, cause locally additional pressure. However, both regions are clearly better 
described by the inclusion of the outer ring. The averaged gas flaring is better described for 
an outer ring mass of 2.1 • 10 10 M Q , while the RC is best described by a twice heavier outer 
ring mass. Given the fact that the outer RC and the gas flaring use data at very different 
longitudes and the northern and southern hemisphere show that the gas flaring is dependent 
on longitude one should not worry about a factor of two. 

3 Summary 

Up to know determinations of the dark matter density typically only considered the rotation 
curve of the inner Galaxy and found values around 0.3-0.4 GeV/cm 3 . New observations of the 
rotation curve challenge this result: the VERA Very Large Baseline Interferometer confirm 
the dip-like structure in the rotation curve around 9 kpc beyond doubt. Such a dip can only 
be caused by ringlike DM structures, since spherical shells do not influence the gravitational 
potential inside. It would be hard to explain the dip by spiral arms or warping of the disk, 
since spiral arms with such high density have not been observed in the outer Galaxy and 
the VLBI measurements used maser light from molecular clouds, which is independent of the 
warping. 

To explain the two broad "dips" at 3 and 9 kpc in the rotation curve two ringlike 
structures are needed. A fit yields a maximum density of the donut-like rings at 4.2 and 12.4 
kpc with a Gaussian width in radius of 2.5 and 3.2 kpc, respectively. The first surprising 
result is not only a reasonable description of the rotation curve including the precise VLBI 
data, but also the fact that the position and size of these rings coincide with the dust ring at 4 
kpc and the Monocerus ring of stars at 13 kpc. The first one is suggestive of a local potential 
well, while the latter is thought to originate from the infall of a dwarf galaxy, as suggested 
by N-body simulations. Of course, it is worth emphasizing that our model of circular ringlike 
structures should not be taken literally, since DM structures will not be perfectly circular. 
However, the data are most affected by the nearby segment of the structure, which seems 
well described by a circular part, at least changing the ellipticity of the rings did not modify 
the results and both, the Monocerus ring and dust ring are largely circular as well. It should 
be noticed that the significant enhancement of the rings in comparison with the smooth halo 
in Fig. 2 may be exaggerated, since recent N-body simulations predict an enhancement of 
the smooth halo in the disk, the so-called dark disks [41], which would reduce the contrast. 
Furthermore, the parametrization with ringlike structures is certainly an oversimplification, 




(2.10) 
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so the error on the local DM density is large and the difference between DM density with 
and without rings is about three standard deviations. 

The second surprising result is the agreement with the s-shaped gas flaring in the outer 
Galaxy, which also needs a ringlike structure with a similar mass at a similar position in 
the outer Galaxy. The "dip" in the gas flaring was determined in such a way that it is 
independent of the warping, basically because the center of the gas layer is not required to 
be at z=0. 

The third surprising result is that a fourfold increase of the local DM density compared 
to a fit without rings is still in excellent agreement with the local surface density, as deter- 
mined from the density distribution and velocity dispersion of the local stars. From these 
measurements it was claimed that our Galaxy needs locally no DM in agreement with the 
decrease of the inner rotation curve, which is suggestive of a potential determined by 
visible matter alone. If the decreasing rotation curve is attributed to the negative contri- 
bution from the outer ring, one obtains a larger local DM density and a correspondingly 
higher surface density. A higher surface density implies a higher vertical force, which only 
changes the slope of the gravitational potential, not its magnitude. However, its magnitude 
is what is probed by the distribution function of local stars. Below z=0.5 kpc the potential is 
dominated by visible matter, so a fourfold increase in the DM density will change the slope 
of the potential smoothly by about 20%, but the data at large z are not precise enough to 
be sensitive to such a moderate change of slope. 

In summary, the higher DM density found in this paper seems inevitably required by 
the structure of the rotation curve using the new VLBI measurements and the gas flaring 
and was found to be perfectly compatible with the surface density. Such a possible increase 
of the DM density in the disk is of prime importance for DM searches, both direct and 
indirect, since the direct detection rate is proportional to the local DM density, while nearby 
ringlike DM structures would change the directional dependence for gamma rays in indirect 
DM searches. 
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